X-Ray Studies of HESS J1809-193 with Suzaku 
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Abstract 



Suzaku observed the region including HESS J1809— 193, one of the TeV uniden- 
Q ■ tified (unID) sources, and confirmed existence of the extended hard X-ray emis- 

sion previously reported by ASCA, as well as hard X-ray emission from the pulsar 
PSR J1809— 1917 in the region. One-dimensional profile of the diffuse emission is 
j- T l'j represented with a Gaussian model with the best-fit a of 7±1 arcmin. The diffuse 

emission extends for at least 21 pc (at the 3a level, assuming the distance of 3.5 kpc), 
^ ! and has a hard spectrum with the photon index of T ~1.7. The hard spectrum 

suggests the pulsar wind nebula (PWN) origin, which is also strengthened by the 
hard X-ray emission from PSR J1809— 1917 itself. Thanks to the low background of 
Suzaku XIS, we were able to investigate spatial variation of the energy spectrum, but 
no systematic spectral change in the extended emission is found. These results imply 
that the X-ray emitting pulasr wind electrons can travel up to 21 pc from the pulsar 
without noticable energy loss via synchrotron emission. 
q\ ', Key words: gamma rays: observations — ISM: individual (HESS J1809— 193) 



stars: pulsars: individual (PSR J1809-1917) — X-rays: ISM 



1. Introduction 

Galactic plane survey with the H.E.S.S. Cherenkov telescope system revealed dozens 
of the new very-high-energy (VHE) 7-ray sources (Aharonian et al. 2005a; Aharonian et al. 
2006a). Many of them have no counterparts in other wave-lengths, thus called "unidentified 
(unID) TeV sources" . Today, about 40 such unidentified TeV sources are known on the Galactic 
plane (Hinton 2007). Most of them are located within a height of ± 1 degree from the Galactic 
plane, and some are intrinsically extended. Despite a large number of intensive studies in the 
last several years, their origin is unclear (Hinton 2007). 

X-ray follow-up observations of the unID TeV sources are now on-going. Although 
supernova remnants (SNRs) or hypernova remnants were suggested to be major counterpart 
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candidates of these TeV unID sources (Yamazaki et al. 2006; Ioka & Meszaros 2009), only a few 
sources have been actually identified as SNRs (Nakamura et al. 2009; Aharonian et al. 2008). 
On the other hand, rather surprisingly, several unID TeV sources have been identified as pulsar 
wind nebulae (PWNe) (Anada et al. 2009; Uchiyama et al. 2009, for example). They seem to 
be rather old, previously unknown PWNe, compared to the PWNe already identified and well 
studied in X-rays. 

The first HESS observations of the region around PSR J1809-1917 were made from 
May through June 2004 as part of the systematic survey of the inner Galaxy (Aharonian 
et al. 2005a; Aharonian et al. 2006a). Because marginal VHE 7-ray signals were detected, 
HESS J1809-193 was observed again in 2004 and 2005, and significant 7-ray emission was 
confirmed (Aharonian et al. 2007). Recent study of this source by HESS was reported by 
Renaud et al. (2008). Fitting the excess map with a 2-D symmetric Gaussian, the best fit 
position and intrinsic source extension (in rms) were determined as (RA, Dec) = (18 h 09 m 52 s , 
-19°23'42") and 0°.25±0°.02, respectively. 

PSR J1809-1917 is a radio pulsar discovered by the Parkes Multibeam Pulsar 
Survey (Morris et al. 2002). The pulsar is located at the position of (RA, Dec) = (18 h 09 m 43fl, 
— 19°17'38") with a pulse period of P = 82.7 ms and the period derivative of P = 2.55 x 
10~ 14 s s -1 . The distance to the source was estimated to be d = 3.5 kpc from the pulsar's 
dispersion measure using the NE2001 Galactic electron-density model (Cordes & Lazio 2002). 
The characteristic age and the spin-down luminosity are rc = 51 kyr and E — 1.8 x 10 36 ergs s _1 , 
respectively. 

The 7-ray spectral analysis performed by Komin et al. (2007) indicated that the spectral 
slope is different between the regions near the pulsar and far from the pulsar. This is the second 
case that spatial variation of the spectral slope is revealed in the VHE 7-ray emission. The 
first case is HESS J1825-137 (Aharonian et al. 2006b), which is largely extended in both VHE 
7-ray and X-ray bands (Uchiyama et al. 2009). 

ASCA observation revealed diffuse, non-thermal emission in the vicinity of PSR J1809- 
1917 (Bamba et al. 2003a). Kargaltsev & Pavlov (2007) detected a bright point X-ray source 
which was positionally consistent with the pulsar PSR J1809— 1917, and resolved the surround- 
ing compact PWN utilizing the very high angular resolution of Chandra. The PWN has a 
"head-tail" profile, consisting of the southern-head, which is coincident to the pulsar, and the 
northern-tail. Kargaltsev & Pavlov (2007) claimed that this cometary morphology is attributed 
to a bow shock created by the pulsar moving supersonically to the southern direction. However, 
it is still unknown whether the faint diffuse emission discovered by ASCA is related to the PWN 
or not. In this paper, we make detailed analysis of the diffuse emission for the first time us- 
ing Suzaku. Suzaku, characterized by the low detector background compared to Chandra and 
XMM-Newton, is very suitable for the analysis of faint and diffuse X-ray emission such as HESS 
1809-193. 
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2. Observations 




Fig. 1. A smoothed excess map measured by HESS around HESS J1809-193 (Aharonian et al. 2007). 
Position of the pulsar PSR J1809-1917 is marked with a cyan cross, whose Chandra image (Kargaltsev & 
Pavlov 2007) is shown in the inset at the top right of the figure. The best fit position and extent of the 
7-ray source is marked with a black cross and a dashed circle, respectively. Two cyan squares indicate the 
two XIS field of views (FOVs) during the Suzaku observations, which are referred to as the North and the 
South pointing, respectively. 

We observed HESS J1809-193 with Suzaku (Mitsuda et al. 2007) in April, 2008. Suzaku 
is equipped with two types of instruments: the X-ray Imaging Spectrometers (XIS: Koyama 
et al. 2007) at the foci of four X-Ray Telescopes (XRT: Serlemitsos et al. 2007) and the Hard 
X-ray Detector (HXD: Takahashi et al. 2007, Kokubun et al. 2007). The observation was 
carried out with two pointings at north and south of the source region (figure 1) in order to 
cover the pulsar and extended VHE 7-ray emission along the direction of the elongated shape 
of the PWN (Kargaltsev & Pavlov 2007). Three XISs (XIS 0, 1, 3) out of four were operated 
in the normal clocking mode with the Spaced-row Charge Injection (SCI) (Nakajima et al. 
2008). We analyzed the data processed by the version 2.2 pipeline. We are interested in the 
spatial variations of a scale of arcminutes. Hence, we concentrated on the analysis of the XIS 
data in this paper, since the HXD does not have a spatial resolution within a field of view of 
~30 arcmin. (FWHM). We applied the standard screening criteria to the XIS data 1 to obtain 
the cleaned event lists. After the data screening, the net exposures was 51.5 ks and 44.2 ks for 
the north and south pointing of XIS, respectively. The Suzaku observation log and exposure are 
summarized in table 1. We used HEADAS version 6.5 software package for the data analysis. 

1 http://www.astro.isas.jaxa.jp/suzaku/process/v2changes/criteria xis.html 
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Table 1. Journal of the Suzaku observations of HESS J1809-193 



North 



South 



Sequence ID 
Start time (UT) 1 
End time (UT) 1 
Aim point R.A. (J2000.0) 
Aim point Dccl. (J2000.0) 
Net exposure (ks) 



503078010 503079010 

2008/03/31 14:06 2008/04/01 16:34 

2008/04/01 16:30 2008/04/02 14:47 
18 h 09 m 37 s .4 18 h 09 m 21 s .0 
-19°21'24" -19°32'02" 
51.5 44.2 



1 Time form of yyyy/mm/dd hh:mm 




Fig. 2. Left: Suzaku XIS image of the region including HESS J1809-193. The pseudo-color represents 
vignetting-corrected, log-scaled intensity levels in the 2.0-10.0 keV band. The yellow cross indicates the 
centroid of HESS J1809-193, and the yellow lines represent 80%, 60%, and 40% contour of the peak of 
the VHE 7-ray (Renaud et al. 2008). Right: normalized 1-dimensional profile of the surface brightness in 
the 2.0-10.0 keV band along the direction from north to south in the green rectangle indicated in figure 2 
with a spatial bin size of l'.O. Solid curve and horizontal dashed line show the best-fit Gaussian profile 
and a constant, respectively. The bright peak corresponding to the pulsar was ignored in this fit (2'. 9 
wide, inner part between the two vertical dashed lines). 

3. Results 

3.1. X-ray Image 

In order to make images of this field, we corrected the vignetting effect by dividing the 
image by the flat sky image simulated using xissim (Ishisaki et al. 2007) after subtracting 
non X-ray background (Tawa et al. 2008). In this simulation, we assumed the input energy 
spectrum as that extracted from the red rectangular region (see §3.2). Hereafter, all images 

4 



are vignetting corrected. Figure 2 shows Suzaku XIS 2-10 keV image of the HESS J1809-193 
region. 

We confirmed largely extended emission reported by ASCA (Bamba et al. 2003a). In 
addition, we found the extended emission has several peaks. The pulsar is the brightest in the 
FOV, and the western edge of the FOV is also significantly bright, which was already hinted 
by ASCA (Bamba et al. 2003a). 

We have determined extension of the diffuse emission as follows: We created a 1- 
dimensional profile of the surface brightness from the rectangle region shown in figure 2 (left) 
which runs from north to south. We selected the direction which enabled us to create the 
longest profile. We did not take a symmetric region around the pulsar, because the faint and 
diffuse emission is significantly asymmetric. Although there is a point source in the south 
of the pulsar at (272.42, —19.43), its flux in the integrated region is only 1 % of the central 
pulsar, thus negligible. The 1-dimensional profile thus created is shown in figure 2 (right). 
Note that the surface brightness is normalized to the peak brightness. We fitted the profile 
with a Gaussian function plus a constant to evaluate its extension. In the fitting, we ignored 
the brightest part around the pulsar with a width of 2'. 9 (corresponding to the point-spread 
function of the XRT). Consequently, we found the diffuse emission extends up to ~ 20' away 
from the pulsar. The Gaussian center was found to be offset by ~3' from the pulsar and the 
rms width to be a = 7' ± 1'. 

3.2. Energy Spectra 

We studied spatial variation of the X-ray energy spectra of the diffuse emission. We 
generated the detector response file and an auxiliary response file using xisrmfgen, and 
xissimarf gen (Ishisaki et al. 2007) and performed model fitting using XSPEC version 12.4.0. 

First we examined energy spectrum from an outskirt of the HESS source region (the red 
rectangular region in figure 2) in order to estimate contribution of the Galactic Ridge X-ray 
Emission (GRXE) on this particular area. The bright source at the south-east corner (circled in 
red) was excluded in this analysis. We subtracted the non X-ray background (NXB) estimated 
using xisnxbgen (detail of the method described in Tawa et al. 2008), and averaged all of the 
available XISs (XIS 0,1,3) data. Because the NXB becomes high above 7.2 keV (in particular 
for XIS1 which has the back-illuminated chip), we used the data only below 7.2 keV. We fitted 
the spectra in 4.0-7.2 keV and 2.0-4.0 keV separately following the reproduction procedure of 
GRXE in Ebisawa et al. (2008). The model adopted is a power-law plus three narrow Gaussian 
(intrinsic width fixed to zero) in 4.0 - 7.2 keV, or four gaussian lines in 2.0 - 4.0 keV. The 
spectra and the best-fit models are shown in figure 3. The equivalent widths of the iron lines are 
summarized in table 2. These equivalent widths are comparable to those of GRXE described in 
Ebisawa et al. (2008). Thus we conclude that the X-ray emission in this region can be regarded 
as pure GRXE and no significant contribution from the pulsar is present. 
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Fig. 3. Left: NXB-subtracted XIS spectrum in the 4.0-7.2 keV band extracted from the region enclosed 
with a red rectangle in figure 2 (left). The best-fit model, an absorbed power-law with three iron lines, 
is also shown in a solid line (each component in a dot-dashed line). The bottom panel shows residuals to 
the best-fit model. Right: the same as the left one, but in the 2.0-7.2 keV band. Narrow Gaussians were 
added to the model below 4 keV. 

Table 2. Center energies and equivalent widths of the iron lines in the background region in comparison of those of the 
GRXE. 





Background rej 


^ion 




Center energy (keV) 


6.44 ±0.07 


6.69 ±0.02 


6.98 ±0.08 


Equivalent width (eV) 


70 (0-140) 


240 (160-550) 


50 (0-170) 


GRXE* 


Center Energy (keV) 


6.41 ±0.02 


6.670 ±0.006 


7.00 ±0.03 


Equivalent width (cV) 


80 (60-100) 


350 (310-390) 


70 (40-100) 



Center energies and equivalent widths determined by the GRXE observation with Suzaku by Ebisawa et al. (2008). 



In the next step, we divided the sky region of the north pointing into a "check pattern" 
as shown in figure 4 (left) in order to find out the possible spatial variations of the spectral 
slope. Here we refer to these regions as the number indicated in figure 4 (left) with the prefix 
"Grid" (Gridl, Grid2, ...). We used a two component model for the fit: an absorbed power-law 
plus the GRXE. We used the model GRXE spectrum (as explained below) to subtract the 
diffuse background. 

The model GRXE spectrum was constructed as follows: The background spectrum was 
fitted in the 2.0-7.2 keV band with the model of an absorbed power-law plus 7 narrow Gaussians 
as explained above. The absorption column density was fixed to 1.0 x 10 22 cm -2 . The best-fit 
parameters, besides the iron line parameters in table 2, are listed in table 3. Normalization for 
each extraction region was adjusted taking account of differences of the vignetting effects, size of 
the extraction region, and exposure time between the north and south pointings. The correction 
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Table 3. Best-fit parameters of an absorbed power-law model with emission lines. 



Model component Parameter Value 



Continuum 


IV tl (W cm 


) 


l.U (nxcaj 




r 




1 40+005 


Emission lines 


Energy (keV) 


Sxv 


2.43 ±0.02 






Sxvi 


2.61 ±0.04 






Ar XVII 


3.13 ±0.02 






Ar XVIII 


o 4 o+0.13 



Note. — Errors represent single-parameter 90% confidence limit. 



factors of the vignetting effect (shown in figure 4 right) were determined by simulation. We 
subtracted the non X-ray background (NXB) estimated by using xisnxbgen. Figure 6 shows 
the XIS spectra (averaged for XIS 0, 1 and 3) and the best-fit models for all the 16 regions in 
the 2.0-10 keV band. Because Grid 1, 4, 13 were illuminated with the calibration sources, the 
data between 5.73-6.67 keV were removed from the fit. Fit results are summarized in table 4. 
The spatial distribution of the spectral indices is shown in figure 5. In order to make statistics 
better, we also fit the combined spectra of Grid 6 and 7 (for the region close to the emission 
peak) and Grid 11, 14, and 15 (for the region far from the emission peak). The results are 
included in Table 4. From these analyses, we concluded that there is no systematic trend of 
the spatial variations in the spectral index in this field of view. 

4. Discussion 

Using Suzaku, we have clearly shown presence of the large scale diffuse emission around 
HESS J1809-193, which was already suggested by ASCA (Bamba et al. 2003a). The extension 
of the diffuse emission is at least ~21' (3cr of the gaussian approximation). The diffuse spectrum 
has the photon index of T ~1.7, which is much harder than those of SNRs with synchrotron 
X-ray emitting shells (r ~2-3), such as SN 1006, RX J1713-3946, and others (Bamba et al. 
2008; Takahashi et al. 2008, for example). Such a hard spectrum, on the other hand, reminds 
us the PWN origin (Kargaltsev & Pavlov 2008). This hypothesis is, in fact, strengthened by 
the very existence of the pulsar PSR J1809— 1917. 

The TeV emission also may have the PWN origin, because it positionally coincides with 
the region including the pulsar and the diffuse X-ray emission, although the center of the TeV 
emission is offset by ~6' (a projected distance of 6 pc at 3.5 kpc) from the pulsar. If both 
the TeV and X-ray emissions are from the same PWN, their origin is accelerated electrons via 
inverse Compton (TeV) and synchrotron mechanisms (X-ray). Typical energies of responsible 
electrons are ~20 TeV for the TeV emission and ~80 TeV for X-rays, respectively. The former 
has much longer synchrotron life timescales than latter (Mattana et al. 2009, for example). If 
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Fig. 4. Left: Suzaku XIS image of the north pointing in the 2.0-10.0 kcV band and designation of the 16 
sub-regions for spectral variation study. Right: Correction factors for each region, which take account of 
the vignetting effects, differences of the exposure and the extraction area between the source region and 
the GRXE region (the red rectangle in figure 2). 




1.2 1.4 1.6 1.8 2 2.2 2.4 2.6 



Fig. 5. Spatial distribution of the spectral indices around PSR J1809-1917. Left: best-fit photon index 
map. Gray scale indicates the photon indices; light gray means the spectrum is hard while dark gray 
means soft. Right: photon indices of each grid with 90% error bars. 
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Table 4. Best-fit parameters of an absorbed power-law model for each grid. 



Grid id 


r 


Flux* (10~ 13 ergs cm -2 s" 1 ) 


X 2 /d-oi. 


1 


1.52 + 0.21 


8.7 + 0.9 


52.2/46 


2 


1.83 + 0.18 


7.8 + 0.6 


70.6/67 


3 


1.79 + 0.21 


5.9 + 0.6 


54.1/62 


4 


1.39 + 0.28 


fi i +o » 


49.1/37 


■5 


1.50 + 0.14 


10.6 + 0.7 


74.0/73 


6t 


2.07 + 0.14 


8.2 + 0.5 


105.0/102 


rt 


1.66 + 0.17 


6.7 + 0.5 


121.7/93 


8 


1 71 +0-40 
1 - 1 i -0.38 


3.0 + 0.5 


59.2/48 


9 


1.71 + 0.18 


7.1 + 0.6 


58.1/61 


10 


1.94 + 0.18 


5.9 + 0.5 


122.9/86 


11 


2 74+0.39 

z - ' -0.35 


2.4 + 0.4 


87.6/70 


12 


i 79+0.38 
1 - ' z -0.37 


2.9 + 0.5 


65.2/47 


13 


1.19 + 0.26 


8.3 + 1.0 


51.4/41 


14 


1.54 + 0.20 


7 9+0.7 
' - z -0.6 


74.3/59 


15 


1.78 + 0.25 


4-4±g : g 


55.5/52 


16 


1.66 + 0.32 


4.9 + 0.7 


60.8/40 


6++7+ 


1.91 + 0.11 


14.7 + 0.9 


235.9/196 


11+14+15 


1.84 + 0.14 


13.7 , '; 


243.2/183 



Note. — Errors represent single-parameter 90% confidence limit. 

* Unabsorbcd flux in the 2.0-10.0 keV band. 

' The emission is dominated by the point source. 



we assume that the electrons which emit TeV gamma-rays have the same age as the pulsar 
(re =51 kyr) and those for the X-rays are very fresh, the positional offset may be explained by 
the proper motion of the pulsar. To explain the offset of ~6', the transverse velocity of the 
pulsar needs to be ~120 km s" 1 , which can be explained with the average transverse velocity 
of radio pulsars (Lyne & Lorimer 1994, ~ 300 km s -1 ). On the other hand, the tails seen 
in the Chandra image (if it is indeed a tail not a jet) suggests that the pulsar is 
moving southward and its velocity vector is at the angle of ~190 200 deg. East of 
North (Kargaltsev & Pavlov 2007). Now if we assume that the pulsar was born 
at the center of the TeV source and moved to its current position, its velocity 
vector would have to be at the angle of 10—20 deg. West of North. Therefore, the 
angle between these two velocity vectors will be ~140— 160 deg., which makes this 
scenario unlikely. Another possible scenario to make such an offset is collision between 
the reverse shock and the PWN (Pavlov et al. 2008) like the case for Vela X (Blondin et al. 
2001; Gaensler et al. 2003). 
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Fig. 6. NXB-subtractcd XIS spectra in the 2.0-10.0 keV band extracted from each grid shown in figure 4. 
Dashed line, dash-dot line and dash-dot histograms indicate the best-fit absorbed power-law model, a 
GRXE model and emission lines, respectively. 

In the young PWNe, it is suggested that the X-ray spectra become softer with 
the distance from the pulsars due to synchrotron energy loss of the accelerated electrons 
(Mori et al. 2004, for example). However, we could not find any hint of such soften- 
ing (Figure 5). This means that electrons do not lose significant energies via the syn- 
chrotron emission. The X-ray size of HESS J1809-193 is ~21' (= 3a), or 21 pc with the 
assumed distance of 3.5 kpc. The synchrotron lifetime of an accelerated electrons (rsyn) is 
rsyn = 6.8kyr( J B/3 fxG)- 3/2 (Esyn/2 keV)" 1/2 , where B and Esyn are magnetic field and mean 
energy of synchrotron emission from accelerated electrons, respectively. In order to explain 
the size of the X-ray diffuse emission, the transport velocity of accelerated electrons should be 
higher than 21 pc/rsyn ~ 3.0 x 10 3 km s" 1 (5/3 fiG) 3/2 (Esyn/2 keV) 1/2 , which seems to be 
very fast for an old PWN, and even comparable to the forward shock velocities of young SNRs. 
It is known that the diffusion coefficient in young PWNe and SNRs is too small to diffuse out 
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to such a large scale, when the magnetic field is turbulent (Shibata et al. 2003; Bamba et al. 
2003b; Bamba et al. 2005, for example). Therefore, the current observation suggests that the 
turbulence of magnetic fields in old PWN systems are smaller than those in young systems. If 
turbulence of the magnetic fields in old PWN systems is smaller than that in young systems, 
such a fast diffusion may be explained, although we have no such a model for old PWNe. In 
supernova remnants, we have some indication that the turbulence becomes smaller when the 
SNR becomes older (Bamba et al. 2005), similar mechanism might work for PWNe. 

We have a similar case, HESS J1825-137, which is another unID TeV source with possible 
PWN origin. Uchiyama et al. (2009) found that the X-ray photon index is nearly constant over 
the emission region, which may suggest HESS J1825-137 also might have small turbulence of 
the magnetic field. Therefore, we suppose fast transportations of accelerated electrons may be 
rather common phenomena in old PWNe. Further similar samples of unID TeV sources may 
confirm this hypothesis. 

5. Summary 

• We observed the TeV PWN candidate HESS J1809-193 with Suzaku, and confirmed an 
extended emission around the pulsar PSR J1809-1917. Size of the X-ray emission is at 
least ~ 21', or 21 pc at 3.5 kpc. 

• The extended emission has very hard nonthermal spectrum with the photon index of 
~1.7. No systematic spatial variation of the photon index is found, which implies that 
accelerated electrons do not lose their energy when they run from the pulsar to the edge of 
the emission region. We have to consider very fast diffusion in an old PWN to reproduce 
such phenomena. If turbulence of the magnetic fields in old PWN systems is smaller than 
that in young systems, such a fast diffusion may be explained. 
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